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THEORETICAL CALCULATION OF THE POWER SPECTRA OF THE ROLLING AND YAWING 

MOMENTS ON A WING IN RANDOM TURBULENCE 1 

By John M. Eggleston and Franklin W. Diederich 


SUMMARY 

The correlation junctions and power spectra of the rolling 
and yawing moments on an airplane wing due to the three 
components oj continuous random turbulence are calculated. 
The rolling moments due to the longitudinal ( horizontal ) and 
normal ( vertical ) components depend on the spanwise distribu- 
tions of instantaneous gust intensity, which are taken into 
account by using the inherent properties oj symmetry oj iso- 
tropic turbulence. The results consist oj expressions for 
correlation functions or spectra of the rolling moment in terms 
of the point correlation functions of the two components oj 
turbulence. 

Specijic numerical calculations are made j or a pair oj correla- 
tion functions given by simple analytic expressions which jit 
available experimental data quite well. Calculations are made 
or jour lift distributions. Comparison is made with the 
results oj previous analyses which assumed random turbulence 
along the flight path and linear variations oj gust velocity across 
the span. 

The rolling moment due to lateral (side) gusts, which is small, 
is expressed in terms of the instantaneous value of the gust near 
the center line oj the fuselage , so that the effect of spanwise 
variation in gust intensity is ignored. The yawing moments 
are considered to be proportional to the rolling moments with the 
constants oj proportionality given by simple aerodynamic 
relations. 

INTRODUCTION 

The gust velocities acting on an airplane flying through 
turbulent air are functions of position or time known only 
in a statistical sense. Consequently, aerodynamic forces and 
moments produced by the lifting surfaces of the airplane 
can be known only in a statistical sense. If the statistical 
characteristics of the turbulence are assumed to be invariant 
with position along the flight path, flight through turbulent 
air may be considered to lie a stationary random process 
and the mathematical techniques developed for such pro- 
cesses (see ref. 1, for instance) may then be used in this 
problem. 

This approach has been adopted in many papers on this 
subject, among them references 2 and 3. Inasmuch as in 
these papers the motions and forces associated with the 
longitudinal degrees of freedom were of primary interest, the 
assumption was made, implicitly, that the gust intensity is 
uniform along the span at any instant. However, for the 


problem of analyzing the motions and forces associated 
with the lateral degrees of freedom, this assumption is 
inadequate, inasmuch as it implies that the vertical and 
horizontal gusts produce zero rolling and yawing moments 
on the wing. This problem has been treated in references 4 
and 5 and elsewhere on the basis of the assumption that at 
any instant the gust intensity varies linearly across the span. 

A fundamental method of accounting for the lift on a wing 
due to random variations of the gust velocities in both the 
flight-path and the spanwise directions is given in reference 6 
for the longitudinal response of an airplane in atmospheric 
turbulence. The approach is based on the assumption that 
the turbulence is axisymmetric (according to ref. 7), so that, 
at any arbitrary time or position in the turbulence, the 
statistical characteristics of the turbulence encountered by an 
airplane do not depend on the heading of the airplane. On 
the basis of this assumption, the variation of gust intensity 
across the span can be related to the variation of the gust 
intensity along the flight path. 

In the present report the approach of reference (i is extended 
to the calculation of the rolling and yawing moments on a 
wing due directly to vertical gusts, longitudinal gusts (herein- 
after referred to as horizontal gusts), and lateral or side 
gusts. These moments are required as a first step in cal- 
culating the motions of a complete airplane in atmospheric 
turbulence; the moments due to the motions caused by these 
input moments can be calculated by conventional methods 
and will not be considered herein. 

In the first part of the report, a theoretical analysis is made 
defining the power spectra of the rolling and yawing moments 
of a wing in terms of the statistical characteristics of the 
atmospheric gust velocities. By using an analytical expres- 
sion to define these characteristics, a numerical solution of 
the lateral moments is presented in the last part of the report . 


a=(3\ll + (kj 2 
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m, m) 
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SYMBOLS 

wing span 
wing chord 

wing mean aerodynamic chord 
complete elliptic integrals of the second and 
first kind, respectively, of modulus k 
longitudinal correlation function for iso- 
tropic turbulence 
Fourier transform of / 


i Supersedes NACA Technical Note 38G4 by John M. Eggleston and Franklin W. Diederich, 1956. 
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lateral correlation function for isotropic 
turbulence 

Fourier transform of g 
indicial-response function of time only 
indicial-response function of time and dis- 
placement 

Fourier transform of two-dimensional corre- 
lation function 


modulus of elliptic integrals, 



reduced frequency, coL/Z7 
modified Bessel functions of the second 
kind 

incomplete modified Bessel functions of the 
second kind 
section lift 

integral scale of turbulence 
rolling moment 
rolling velocity 
dynamic pressure 

yawing velocity (used only in stability 
derivatives) ; linear displacement between 
any two points 


wing area 
time 

mean forward velocity 
displacement along the flight path 
component of airplane velocity along posi- 
tive } r -axis 

three components of gust velocity (see fig. 


1(a)) 

reference axes (see fig. 1(a)) 
chord wise distance 
span wise distance 


nondimensional spanwise coordinate, 
angle of attack, radians 


JL 
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span influence function 
integral weighting function 


dummy variable of integration, y 2 *—yi* 
atmospheric density 
dummy variable of time 


circular frequency, 27 r/Period 

. Rolling moment 
rolling-moment coefficient, 

„ . Yawing moment 
yawing-moment coefficient , qSb~ 
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Subscripts: 

o 

g 

A bar over 


correlation function 
power spectral density 

trim value 
gust component 

t quantity denotes the mean value of 


the 


quantity. The absolute value of a quantity is denoted by | 


THEORETICAL ANALYSIS 
PRELIMINARY CONSIDERATIONS 

In this section expressions are derived for the power 
spectra of the rolling and yawing moments of an unswept 
airplane wing or thin lifting surface of arbitrary plan form 
due to flight through random atmospheric turbulence. 
Essentially, the procedure consists of expressing the rolling 
moment at any arbitrary position along the flight path in 
terms of the gust velocity at that position, establishing the 
correlation function between the rolling moments at any 
two points along the flight path, and transforming this cor- 
relation function into an expression for the power spectral 
density. The power spectrum of the yawing moment is then 
related to that of the rolling moment through simple aero- 
dynamic relationships. 

Assumptions. The following assumptions are made in 
the analysis: 

(1) The turbulence is homogeneous and isotropic; that is, 
the statistical characteristics of the turbulence are invariant 
under a translation or rotation of the space axes (although 
the results obtained for the vertical component of turbulence 
require only the somewhat less restricting assumption of 
axisymmetry). 

(2) Time correlations are equivalent to space correlations 
along the flight path — an assumption usually referred to as 
Taylor’s hypothesis. (See ref. 7.) 

(3) The chordwise penetration factor (the indicial- 
response influence function) for the rolling and yawing 
moments can be expressed as a product of a function of 
distance along the flight path (or time) only and distance 
along the span only. 

(4) The wing considered herein is relatively rigid and, as 
a result of the turbulent velocities, performs small motions 
about a mean steady flight condition. 
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The implication of these assumptions and the limitations 
they impose on the results of the analysis are discussed in a 
subsequent section of the report. 

Coordinate system and gust components. — The system of 
axes and the local velocity field relative to the lifting surface 
are shown in figure 1 (a). The velocity at each point in the 
field is resolved into components lying in the three planes of 
an orthogonal set of axes, the X-axis of which is tangent at 
every point to the flight path. Throughout this report these 
three components are designated as follows: The component 
alined with the X-axis is referred to as the horizontal gust 
u g ; the component alined with the F-axis is referred to as 
the side gust v g ; and the component alined with the Z-axis 
is referred to as the vertical gust w g . 

As the wing moves through the local velocity field, the 
random variations in the horizontal and vertical gust com- 
ponents are defined both in the flight-path direction and in 
the spanwise direction at every position along the flight path. 
Random variations of these gust components across the 
chord are taken into account by indicial-response functions 
and, hence, need not be considered separately. 

The side gust component of the gust velocity field is 
treated in only a limited manner. Neither the chordwise 
nor the spanwise variations of v g are considered along the 
flight path; rather, v g is assumed to act on the wing as a 
point velocity with a variation only along the flight-path 
direction. Contemporary aircraft exhibit such wide varia- 
tions in distribution of dihedral across the span that it is 
doubtful that a generalized analysis could be utilized. The 
point or centroid analysis should be fairly accurate when the 
dihedral distribution is predominant over only a small section 
of the span near the fuselage. Such a distribution is ex- 
hibited by an unswept wing with zero geometric dihedral 
mounted very high or low on a fuselage. For a wing with 
zero aerodynamic dihedral, this component could be neglected 
completely. 

Definition of gust correlation functions. -In order to define 
random variations of the gust velocities both along the 
flight path and across the span of the wing as it moves 
through the turbulence, it is necessary to define the correla- 
tion between any two velocities in the gust field through 
which the wing passes. The space correlation function of a 
velocity u is defined in terms of the distance r as 

1 r x 

^ u (r)=lim 2XJ (!) 

Von Karman and Howarth (ref. 8) have shown that, in 
homogeneous isotropic turbulence, the correlation between 
two velocity vectors a distance r apart can be defined in 
terms of two scalar functions /(r) and g(r) and that this 
relationship is invariant with respect to rotation and reflec- 
tion of the coordinate axes. These one-dimensional corre- 
lation functions relate the paired velocity components 
obtained by resolving the velocity vector at any two points 
a distance r apart into two parts: The pair lying along the 
straight-line path between the points are known as the 
longitudinal components and the pair normal to the straight- 
line path are known as the lateral components. These two 


pairs of components are pictorially shown in figure 1(b)- 
Such velocity components may be measured in wind tunnels 
downstream of a grid mesh. (See ref. 9.) 

In reference 8, it is further shown that these correlation 
functions are interrelated by the differential equation 

5^T+f(r)=S(r) ( 2 ) 

By defining the variable r in the coordinate system of this 
report and using the correlation tensor of reference 8, a two- 
dimensional analysis of the turbulence as it affects the wing 
may be made in terms of f(r) and g(r). The variable in the 
correlation functions of the horizontal and vertical gust com- 
ponents in the two-dimensional XF-plane of the wing is 
given simply by 

r=V(Az) 2 + (A?/) 2 = V(Cr) 2 + (Ay)* (3) 

The correlation function of the horizontal gust components, 
as derived from the correlation tensor of reference 8, is 
defined in t'prms of x- and ^-components of the present 
analysis by the formula 

(iiw 1 } <4) 

The relationship between the components is shown schemati- 
cally in figure 2 (a). 

In a like manner, the correlation function of the vertical 
gust components affecting the wing, given in t erm s of the 
mean-square value of the vertical gust velocity w 2 , may be 



• — - • — * Longitudinal components, f (r) 

H r * — -I 

(b) Lateral components, g (r) 

(a) Wing passing through three-dimensional turbulence. 

(b) Components of turbulence as a function of distance r. 

Figure 1 . — Sign convention and stability axes of a wing passing 
through a turbulent velocity field. Arrows denote positive direction, 
where applicable. 
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seen to be simply 

y Wg (Ax,Ay)=w g 2 g [V (Ax)' 2 -\- (A//)-J (5) 

For the case of side gusts acting on a wing, the correlation 
function would be defined in terms of Ax and Ay if the span- 
wise correlation were considered. (See fig. 2(h).) Inasmuch 
as the side gust is considered to act only at a point on the 
span, Ay is zero, and the correlation function for the side 
gust in terms of its mean-square value becomes 

V Vg {Ax)=v?g(Ax) (6) 

Although the mean-square value of each of the three gust 
components is given separate identity, under the assumption 
of isotropy 

% 2 =v?=w? (7) 

With l lie gust-velocity correlation functions thus defined, 
the forces and moments due to antisymmetric components 




(0) 




(a) Horizontal gust components. 

(1)) Side gust components. 

Figure 2. — Schematic drawing of the relationship between the compo- 
nents of horizontal and side gusts at any two arbitrary points. 


of the gust-velocity field acting on a wing passing through 
that field may be derived in terms of these correlation 
functions. 

ROLLING MOMENT DUE TO GUSTS 

Vertical gusts. The instantaneous wing rolling moment 
due to vertical gusts can he written in terms of an indicial- 
response influence function h'(t,y) as 

J 00 /* 6/ 2 

h'(t u y)w g (t — f u y)dydt , (8) 

- oo J -6/2 

According to assumption (3) of the section entitled “Pre- 
liminary Considerations” (see also the argument presented 
in ref. 6), the function h'(t,y) can be expressed in the form 

h'(f 1 ,y)=h(t 1 )y(y) (9) 

where y(y) is a steady-state span influence function and h(t\) 
contains the unsteady-lift effects. The rolling moment can 
then also be written as 

/»00 /*6/2 

M x (t)= h{ti)dti y( y)w g {t—t 1 ,y)dy (10) 

J - a> J — 6/2 

If the correlation function for the rolling moment is calcu- 
lated from equation (10) and a power spectrum for the 
rolling moment is obtained by taking the Fourier transform 
of the correlation function, the resulting expression may be 
shown to consist of a product of two functions: One function 
is the result obtained from quasi-steady considerations alone, 
and the other is the absolute squared value of the unsteady- 
lift function for sinusoidal gust penetration such as that 
given by Sears in reference 10. Consequently, considera- 
tion will he confined to an analysis using quasi-steady ex- 
pressions for the rolling moment ; that is, the lag in buildup 
of lift across the chord of the wing due to the gusts is not 
included. 

In quasi-steady flow, the rolling moment of a wing due 
to a variable angle-of-attack distribution across the span is 
given by 

M x =qSbC l 

C&/2 

= [l(y)] a=a * w ydy (li) 

«/ — 6/2 

where section lift 

l(y)=ci(y)qc(y ) 
and local angle of attack due to gusts 

a g = wJU 

Now, one theorem of linearized airfoil theory states that 
the lift (or rolling moment) on a wing due to an arbitrary 
spanwise angle-of-attack distribution is equal to the integral 
over the entire wing of the product of the spanwise lift dis- 
tribution due to a unit constant (or linearly varying) angle 
of attack and the given arbitrary angle-of-attack distribu- 
tion. Hence, the rolling moment is also given by 

J 6/2 

[l(y)] a=v a g (y)dy (12) 

- 6/2 
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This theorem is valid not only in steady but also in indicia! 
flow. (See reciprocity theorems of ref. 11.) 

When the indicated substitutions are made, the rolling- 
moment coefficient along the flight path is 


r M G* I ' 11 ' 2 f ct(y)c(y) T~" w g (a;,y) 

Cl{x) ~ b> U 


dy 



X 


y(y*)w g (x,y*) dy* 


(13) 


where y*= ^ and the steady-state lift distribution 


J’ 


(14) 


pertains to a linear antisymmetric angle of attack across the 
span. It may be seen that, by virtue of its definition, 
7 ( 7 /*) must satisfy the relation 

j y(y*)y*dy*=2 (15) 

Horizontal gusts. -In analogy to the analysis of the pre- 
ceding section, consideration will be confined to the quasi- 
stcady case. When stability axes are used, a change in 
forward velocity at any spanwise station increases the mag- 
nitude but does not change the direction of the lift and drag 
vectors. Thus, the horizontal-gust contribution to the 
dynamic pressure is 

H{y)=\ P {[u g {y)+UY-U>} 

=^p(u g 2 +2Uu g ) 



under the assumption that u g <^U. When this linearized 
approximation is used, the lift on each section is proportional 


to the local angle of attack: 

(W =9a 

l(w g ) " °w g 

The rolling-moment coefficient due to horizontal-gust ve- 
locities is thereby defined as 

C, (x)= 2 J - j ^ y{y*)u g (x,y*) dy* (16) 

where now 



The only difference in evaluating 7 ( 7 /*) for horizontal and 
vertical gusts lies in the definition of the parameter having 
a variation of y* across the span; for the vertical gust, that 
parameter is taken as the additional angle of attack and, for 

Us 

the horizontal gust, that parameter is 2a v p* The condition 
that 



remains unchanged. 

ANTISYMMETRIC SPAN INFLUENCE FUNCTION 7 (y*) 

The antisymmetric span influence function 7 ( 7 /*) is defined 
over the span so that any given distribution of 7 ( 7 /*) will 
produce a unit rolling moment. These distributions refer to 
the span loading due to a linear angle of attack a=y* for 

the vertical gust or a linear leading-edge velocity 2a 0 jj = V* 

for the horizontal gust. Four basic variations of 7 ( 7 /*) 
have been considered with the proper constants so that 
equation (15) is satisfied. The equations for the 7 ( 7 /*) vari- 
ations considered are given in table I and plots of these 
variations are shown in figure 3(a). The names given to the 
four distributions obtained by rolling the wing refer to the 
distributions which would be produced by a uniform angle 
of attack. 


CORRELATION FUNCTION OF THE ROLLING MOMENT 


Vertical gusts. — The autocorrelation function of 
path of the wing is defined as 


the rolling moments due to vertical gusts at any two stations along the 


^cix 2 — £i) = lim 

1 AT->oo 


X-XX 


C,(x 2 )C t (xi)dx { 


(18) 


With the substitution -of the expressions for Ci(x 2 ) and C t (x i) as given by equation (13), the correlation function of equation 
(18) becomes 

Zi)=hm j ^ : iqzp J i 'Y(Vi*)™*(xuyi*)dy 1 * j y(y 2 *) w g (x 2> y 2 *) dy 2 * dx y 

= I 6C 2 / j 7 ^*) 7 ^ 2 *) hm _ x w s( x uy*) u>g(x 2 ,y 2 *) cfoi dyf dy 2 * 

Ci* n n 

= idp J _! J i 7 (yi*) 7 ( 2 / 2 *) '*w g (x 2 -x l ,y 2 *-y l *) dyfdy 2 * 


( 19 ) 
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where it is assumed that the functions are convergent under 
either order of integration. An expression is thus obtained 
for the correlation function of rolling moment in terms of the 
correlation function of vertical-gust velocity. In equation 
(19), 

i r x 

^tv g (x 2 —x l ,y^—y l *)=hm — ^ w g {x u y^)w g {x 2) y 2 *)dx x 

V->oo J -X 

is the same as the two-dimensional correlation func- 
tion defined earlier as equation (5) with x 2 — x x =Ax and 

y °- yi - Ay ~bj2 

By the proper substitution of variables, the double inte- 
grals of equation (19) may be separated into the single 
integral of the product of the integrated weighting functions 
of 7 (y*) and the correlation function ^ Thus, with the 
substitution of 

x 2 — X\ = Ut 


equation (19) becomes 


V2*~yi*=V 


where 


and 


VcffJr) j r (v) 'i'w/U T,v) dv 

r«-j: V y(yi*)y(yi*+v)dy l * 

{V ( 


*..(Vr,n) =< g | A / (UrY + 


( 20 ) 


( 21 ) 


( 22 ) 


Equation ( 22 ) may be recognized as being equivalent to 
equation (5). 

Horizontal gusts. -In an identical manner, the autocorre- 
lation function of rolling coefficient due to horizontal gusts 
at any two stations along the path of the wing is derived by 
use of equations (16) and (18): 

1 C x 

*cXz 2 — x\)= hm tt-t? Ci(x 2 ) Ci{x x ) dx x 
x_>co 2AJ-X 

<*o 2 c t 2 n r 

= 7 ^*) ^u g {x 2 -x u y 2 *-y^)dy^dy 2 * 

With the same change of variables as in the preceding section, 


{Ur) = °2[/2 ~J () r 00 *u g (Ur,v) dy 


(23) 


where 


^uXUr,y)=U 2 h 




(UrY + (^) 

m 

(UrY + (^j 


^ 4 w*[V^+(W] 


and equation (24) is now the equivalent of equation (4). 
The integral weighting function T{rj) is the same for both 
the horizontal- and the vertical-gust contributions to their 
rolling-moment correlation functions. 

INTEGRAL WEIGHTING FUNCTION T( v ) 

The integral weighting function T{y) as defined by equa- 
tion ( 21 ) has been evaluated for the four distributions of 
7 (y*) given in table I. These values are listed in table II 
and plotted against y in figure 3(b). It may be shown that 
the nature of the function is such that the relationship 


j: 


T{y)dy=0 


(25) 


must be satisfied for any variation of T which pertains to an 
antisymmetric variation of 7 (?/*). In table II the elliptic 
distribution is given in terms of K(k) and E{k), which are 
complete elliptic integrals of the first and second kind, 

9 ~ 

respectively, of modulus k=~ — • The derivation of the 

2 + 77 

elliptic weighting function is included in the appendix of 
the report. 

POWER SPECTRA OF THE ROLLING MOMENT 

The power spectrum of the rolling-moment coefficient C t 
is defined as the Fourier transform of the autocorrelation 

TABLE I 

VARIATION OF y (y*) 


Distribution 

7 (y*) 

Rectangular 

6 y* 

Elliptic 

~ ?/*V 1 —y* 2 


7 r 

Parabolic 

loy*(l-y* 2 ) 

Triangular 

24y*(l-|y*|) 


TABLE II 

INTEGRAL WEIGHTING FUNCTION rfo) 


Distribution 

T(v) 

Limits 

Rectangular 

6(4 — 677 + ?7 3 ) 

0=, = 2 

Elliptic 

512 

^(2 + ,) [4,(^-3,- < l)IC(fc) + 
(4+V -,«)£(*)] 

0S,^2 

Parabolic 

'4 (64 - 336, 2 + 280, 3 - 42,5 + 3,7) 

Zo 

0^,S2 

Triangular 

OQC 

^ (2 - 1 0,2 + 5,2 + 5, 4 - 3,5) 

15 

02=, = 1 


OOO 

^(8 - 20, + 1 0,2 + 5,2 - 5, 4 + ,5) 
15 

1 =, = 2 


(24) 
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(a) Variation of span (b) Variation of integral 

influence function. weighting function. 

Figure 3. — Variation of span influence function y and integral weight- 
ing function r for four types of loading. 


function of Ci : 




» * Ci {Ut) e ~‘Tf Ur d(U t) 


(26) 


Thus, for two of the three components of turbulent gust 
velocities, the power spectrum of the rolling-moment co- 
efficient is dependent on the integration of a function of the 
lifting distribution of the wing times a function which repre- 
sents the Fourier transform of the correlation function of 
the vertical and horizontal gust components over the wing 
span. 

As previously stated, these results are based on quasi- 
steady considerations. Unsteady-lift effects can be taken 
into account simply by multiplying the power spectral 
density of the rolling moment due to each gust component 

by the function |^^)| ' where v is the Sears function 
given in reference 10. 

APPROXIMATION FOR SIDE GUSTS 

As pointed out previously, the side gust is treated here 
only in an approximate manner; that is, the spanwise effect 
is neglected. Based on this approximation, the rolling- 
moment coefficient is defined as 


0^=0,/-^- 


(31) 


For the vertical gusts, the power spectrum of the rolling- 
moment coefficient may be found by substituting the 
derived relationship for Vc^Ur) given by equation (20) into 
equation (26): 


The correlation function is defined by 

*c ( (Ur)= 6 ^' V g(Ur) (32) 


g^ipj e U j ^(v)^u, g (UT, y)dyd(Ur) 


and the power spectrum is defined by 


C t 2 (• 2 -i^-Ur 

= 8irf/\j„ r ^J 6 ' v)d{Ur)dy 

= 87rU 3 jo r(v> Iw (jT V ) dv (27 ) 

Changing the order of integration here is permissible in- 
asmuch as the integrals of the correlation function of w g are 
convergent in both Ur and y. The integral I w is defined as 

u’ vj- e ~‘ v T ’*v> g ( Ur ’ r i) d( PT) ( 28 ) 

Similarly, the power spectrum of rolling-moment coefficient 
due to the horizontal component of gust is obtained from 
the substitution of equation (23) for the term ^c ; (Z7t) 
appearing in equation (26) : 

* c ‘ ( “ )== l Iu (U’ v ) dv ( 29 ) 

where the integral I u in equation (29) is defined as 

h (jj> v)=f_ S* y) d(U r) (30) 


< I’c 1 (co)=%^ G(u) (33) 

where 

t7(co)=-Lj“ e - i! b Ut s (Ut) d(U t ) (34) 

RELATIONS BETWEEN THE YAWING AND THE ROLLING MOMENTS 

No attempt is made herein to calculate directly the yawing 
moment due to atmospheric turbulence. Because of the 
more complicated nature of the phenomena which give rise 
to drag, as compared with those which give rise to lift, such 
an undertaking would be quite difficult. Furthermore, in 
view of the fact that the yawing moments on the wing due 
to turbulence are relatively small, a detailed analysis would 
not generally be warranted. In this section, therefore, an 
approximate procedure is outlined for obtaining the yawing 
moments from the rolling moments. 

The yaw'ing-moment coefficient due to sideslip can be 
expressed in the form 

C n =C nfi (a)P 

where, in this case, a is the sum of the trim angle a 0 and the 

yj 

instantaneous mean vertical-gust angle jjj and where /3 is 


435574—58 2 
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the instantaneous mean 


side-gust angle 


w 


so that 



V d« ) a U U 


where the second term is of higher order and is neglected. 
Similarly, differences in v R along the span give rise to higher 
order terms. 

The rolling moment can be expressed in the same form, so 
that the relationship between the yawing and rolling mo- 
ments due to side gusts is given by 



Actually, this contribution to the yawing moment is generally 
negligible and is included here primarily for the sake of 
completeness. 

For the yawing moments due to vertical and horizontal 
gusts, similar reasoning may be employed. The yawing 
moment in these eases arises from the antisymmetric part of 
the instantaneous angle-of-attack distribution due to turbu- 
lence, as does the rolling moment, so that the two moments 
may be expected to be approximately proportional to each 
other; that is, 



In essence, these relations imply that the yawing moment 
due to a given instantaneous spanwise gust distribution is 
the same as the yawing moment due to a linear gust distribu- 
tion which has the same rolling moment. The deviation of 
the actual distribution from a linear one results in small 
differences in the vortex field and, thus, in small differences 
in the induced downwash. These differences lead to a con- 
tribution to the yawing moment which is believed to be small 
and, hence, has been ignored. 

In terms of their power spectra, the yawing moments are 
defined as 


**.(“) 

= 1 

(&) 

. 2 ,A 

4 V,(«) 


1 

\u K 






(&) 

1 4v ( M 

► 


v g 

r'w 



4>c„h| 

-1 

m 

1 $c,(") 



"d 

(A 'pJ 

“o 



The power spectra of the rolling moments are defined in the 
preceding sections. 


APPLICATION 

APPROXIMATIONS TO THE ONE-DIMENSIONAL (POINT) CORRELATION 

FUNCTIONS 

In order to evaluate the effects considered in the preceding 
part of the report, calculations will now be made by using 
the results derived therein. These calculations will be based 
on a simple analytical expression for the longitudinal point 
correlation function which lias been suggested in reference 
12 on the basis of measurements in wind tunnels: 

f(r)=e L (39) 


where L is the longitudinal scale of turbulence defined for 
any longitudinal correlation function /(r) by 



The characteristics of clear-air turbulence measured in 
the atmosphere (ref. 13) may be shown to be reasonably well 
represented bv equation (39), with a value of L of approxi- 
mately 1,000 to 2,000 feet. There are some theoretical 
objections to this function primarily the fact that it has a 
nonvanishing slope as r — >0 and, hence, that the associated 
power spectrum does not decrease rapidly enough for very 
short wavelengths. These conditions imply that the mean 
square of the derivative of the gust velocity with respect to 
the space coordinate is infinite. However, from available 
measurements on atmospheric turbulence, it appears that 
equation (39) remains valid to distances which are small 
compared with the span of the airplane (on the order of 
several inches), and the behavior of the spectrum at very 
short wavelengths is relatively unimportant because airplanes 
cannot respond to them to any appreciable extent. There- 
fore, in the absence of more reliable information all calcula- 
tions described in this report are based on equation (39). 

The corresponding lateral correlation function related to 
f(r) by equation (2) is found to be 

(4i) 

A plot of the functions given by equations (39) and (41) 
is shown in figure 4. Their respective power spectra, denoted 
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Figure 4. — Variation of lateral and longitudinal correlation functions 

with \r\/L. 


by G(Jc') and F(k') where k '=“{-, are given by 


F(lc') = 


2 L 1 
■kU i I (k-'y- 


(42) 


L l-)~3(i / ) 2 


G(k')=^r ^ 


TU[i+(k')r 


(43) 


These power spectra are plotted to a logarithmic scale in 
figure 5, where it may be noted that the asymptotic slope as 
>oo has a value of —2.0. 



Figure 5. — Power spectra of lateral and longitudinal components of 
isotropic atmospheric turbulence. 


CALCULATIONS FOR VERTICAL GUSTS 

Rolling-moment correlation function. — When equation (41) is substituted into equation (22) with 
the correlation function defined by equation (22) becomes 



(Ur) 2 + 



T ,Tj) =W g 2 |j — .yj 


kVxWh- 


m 


-zV<tw+(£)’ 


(44) 


Inasmuch as the evaluation ol the rolling-moment correlation function, as such, is not necessary to the analysis of this report, 
only limited consideration is given to the calculation of autocorrelation functions. Plquation (20) lias been evaluated in 
closed form for the case of the rectangular distribution of the span influence function 7 ( 7 /*) as given in tables I and II: 




dr) 


3w g 2 Ci 2 


y 4(X 2 +6+j8' 2 ) Vi XV + ( 1 2X 2 - 0 /2 X 2 +24 + 1 ) c 

(4X 3 + 1 2X 2 — 3/3 ,2 X 2 -|- 24X + 24) e _x -f-/3' 3 X [2?i (/S' ,X) — X/? 0 (/3' , X) ] ^- 


— V/9'2+X 2 


- X 2 


( 45 ) 
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where Ko and K x are defined in reference 6 as incomplete modified Bessel functions where 


and 


K 


( / 3 ',x ) =J o 


sinh-'^- 


-XCO8U0 cosh vQde 



x= 


Ur 

L 


(46) 


These two parameters represent the ratios of the distances b and Ut to the integral scale of turbulence L. The parameter 0' 
reflects the size of the wing span relative to the characteristic size of the turbulence and, as such, is one of the more important 
parameters appearing in all the calculations involving spanwise correlation. It effectively scales the magnitude and shape 
of the correlation functions and power spectra and, in the limit as 0'— >0, the equations for the antisymmetric moments like- 
wise go to zero inasmuch as no rolling or yawing moment will exist when a finite span shrinks to a point. 

The parameter X is a measure of the flight-path distance relative to the characteristic size of the turbulence and, in the 
limit as X— >0, the correlation function must reduce to the mean square value of the rolling-moment coefficient; hence, 


Oz 2 =*c,(X= 0) 

3 ^}C X 2 

- g/4 jj r m /3 + 12^ 2 +24^+24)^-^+g /3 -24] 


(47) 


Inasmuch as no adequate tables appear to be available for the functions K 0 and Ki, a numerical evaluation of equation 
(45) has not been made. However, an analysis of this correlation function with other approximations for j(r) indicates that 
the effect of span loading is minor and that a reduction in /3' attenuates the correlation function. 

Evaluation of I v (^>77^. — For the vertical gust component, the integral definition of I w (^^ j7 l^ is given by equation (28) 

whereas '&w g (UT,r]) is now defined by equation (44). The indicated integration may be performed in closed form as a function 
of rj and the reduced frequency parameter k' . Thus, 


7 ”(^) =W ‘ 2 I1 d(Ur) 

=1P lJ-m rfi'n 1 ^[!+3(A:0 2 ] rff 
_ e [^ Vl+(fc/)2 ]+ -i r +(F)T- [t^ 8 ] 


(48) 


where 


V=— 
A. jj 


and Ko and K x are modified Bessel functions of the second kind of order 0 and 1, respectively. 

A plot of equation (48) is shown in figure 6 as a function of the frequency parameter k ' , for a range of values of P'ti/2 
from 0 to 1.0. Although the physical significance of the function I w is rather obscure, the plots are useful in the subsequent 
numerical integration of the product of I w and T. 

Power spectrum of rolling moment. — In general, the analytical solution of equation (27) for the power spectrum of the 
rolling-moment coefficient due to vertical gusts, when possible, is a tedious process. Numerical integration by means of 
either Simpson’s rule or some integration process of higher order is generally preferable to integration in closed form. How- 
ever, the analytical evaluation of equation (27) for the case of a wing with rectangular span loading is given here in order to 
illustrate some of the characteristics of the equations. After the indicated substitutions are made, equation (27) becomes 




^(*0 = 


bwyLCi 


8tU 3 


* n 


(j v + v 3 ) 


r-m 

Li+m 2 


K o- 


^[1+3(A:') 21 
[l + (fc') 2 P 


Ki 


dr] 
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for which the integrated solution is 


\&w*LCi 2 / ra 

*c, (F) = 7rC7 3 ffi 4 [1 + ( ; 0T (« 3 ^') 2 J n K 0 (x)dx+ { a<+ 1 6a 2 [1 - (F) 2 ] } K 0 (a) + 

{ 2a 3 [3- (F) 2 ] + 32a [1 - (F) 2 ] }#,(«) +2a 2 [1- 3(F) 2 ] -32 [1 - (F) 2 ]) (49) 


where -|-(F) 2 , F=^> and K 0 (a) and K, (a) are 

modified Bessel functions of the second kind 2 of argument 
a. Equation (49) is plotted in figure 7 (a) as a function of 
k' for a range of ft ' between 0.03125 and 1.0. 

For small values of frequency co (and hence k') or scale 
factor p', equation (49) becomes poorly behaved because 
the solution takes the form of small differences of high-order 

2 Values for the integral of Ko may be found in several publications, one of which is reference 
14, table 2 (Zahlentafel 2). A comprehensive listing of other available mathematical tables 
including these Bessel functions is given in reference 15. 


terms. The reason for this may be seen by expanding the 
Bessel functions in their power-series form and grouping 
like powers of the variable a. The coefficients of the first 
three terms of the power series a -4 , a -2 , and a 0 (which are the 
predominant terms for values of a<Cl) are identically zero. 
Under these conditions, small computing errors or the lack 
of significant figures will cause large inaccuracies in the 
numercial evaluation of the function. 

The difficulties just described may be overcome somewhat 




k‘ 


(a) Rectangular span loading. 

Figure 7. — Power spectra of rolling moment of wing due to vertical 

gusts. 


a number of values of 0V2. 
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by evaluating equation (49) for the limiting case of &' = (): 



( 50 ) 


When the Bessel functions are again expanded in powers of 0', only several terms are needed to evaluate the function at 
small values of 0'. As before, the coefficients of all negative orders and the zero order of 0' are identically zero. 

The physical necessity that, as the span b approaches zero, the expression for the power spectrum of the rolling-moment 
coefficient must also approach zero is satisfied by equation (49) inasmuch as the lowest order term with a nonzero coefficient 
appearing in the equation is a 2 (as pointed out above) ; that is, for 6— >0, 


In order to compute $c l for the other three types of distribution of wing loading given in tables I and II, a numerical- 
integration process involving Simpson’s three-point rule of integration was employed. The power spectra thus obtained 
are plotted in figures 7 (b), (c), and (d). This method was also used for the rectangular lift distribution arid was found to 
give good agreement with the analytical results. 

It is of interest to note that whereas the power spectra of the vertical gust approach a logarithmic decrement of —2 
(see fig. 5), the rolling-moment power spectra shown in figure 7 approach a decrement of —3. At the low-frequency end of the 
spectrum (long wavelengths) the power appears to approach a constant which is zero only when 0', the ratio of span to scale 
of turbulence, is zero. 

Some simplified approaches to the calculation of the rolling power of gusts (for example, ref. 4) lead to the result that 
the spectrum of the rolling power of the vertical gust appears as the first derivative (slope) of the vertical-gust spectrum. 
As may be seen from figure 7, such an approximation is justified only in a very small band of frequencies for wings having 
small values of 0' . 


Rolling-moment correlation function. — When the expressions for f{r) and g(f) given by equations (39) and (41) are sub- 


The correlation function of rolling moment is obtained by inserting equation (51) into equation (23) and integrating. For 
a rectangular distribution of 7 (y*) } 


4> Ci (o>) ~ (Constant) a 2 =0 


CALCULATIONS FOR HORIZONTAL GUSTS 


stituted into equation (24) with r=^j (Ut) 2 +(^^ , the one-dimensional correlation function for horizontal gusts becomes 



(51) 





(52) 


and 



( 53 ) 
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As in the case of the vertical gust, the correlation function has not been calculated for the other three distributions of y (y*) 
for the reasons already given. 

Evaluation of I u (jj> — The evaluation of I U} as defined by equation (30) for the case of the horizontal gust, is 

given by the expression 

1 1 rW (t7r)2+ © 2 


7 " (jf 


-i-Ur 

e u 


1— 


2 L 




d(Ur) 


r /& 

9 


=Uu 2 L -< 


K ‘ [§ Vi + w 1 5 ]-(&)’ k. [§ vi + W) 


'J 


(54) 


(joL 


P'v 


where and K 0 and K x are the modified Bessel functions of the second kind of argument ~ ^ l + (&') 2 . 




k' 


(b) Elliptic span loading. 
Figure 7. — Continued. 


(c) Parabolic span loading. 
Figure 7. — Continued. 
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(d) Triangular span loading. 
Figure 7. — Concluded. 


The function given by equation (54) is plotted against k' 
in figure 8 for values of $’r\l2. Values taken from this plot 
may be used in the numerical integration of the power 
spectra of the rolling-moment coefficient. 

Power spectrum of rolling moment. — The power spectrum 
of the rolling moment due to the horizontal components of 
turbulence acting on the wing has been determined by using 
the expression for I u obtained in the preceding section and 
the four distributions of the parameter T(rj) given in table II. 
The integral of equation (29) has been evaluated numerically 
for all four cases of load distribution, and the resulting 
variations of the power spectrum with frequency and /3' are 
plotted in figure 9 . In addition, the analytical solutions for 
the cases of rectangular and parabolic distributions are given 
here and their numerical values were checked against those 
obtained by the numerical-integration process. By use of 
equation ( 54 ), the solution for the rectangular case is found 
to be 




ZUgljX-oCi* 


T 2 2 CV0 K 

Jo (4 - 6 ^uvt+(FF 1 \ 2 / . 


drj 


2Au g 2 LoL 0 2 C x 2 


r 

"ir[/ 3 a 4 U + (&') 2 ] L Jn 

16a(a 2 +4) ^(a)— 64^| 


K 0 (x) dx+a 2 (3a 2 +32)K 0 (a) + 



for a range of values of p'v)/2. 


where a=/3'VT+(F) 2 and The analytical solution 

for the parabolic distribution is given by 
^ I 2 ^ 2 

4>c ; (&') = 7 ^ 73 ^ ^ 2 ] [ (209a 6 + 25 ,344a 4 +331,77 6a 2 ) K 0 (a) + 
(a 7 + 3, 424a 5 + 133,632a 3 + 663,552a) K t (a) + 

(a 7 +63a 5 ) P K 0 (x) <fa-l,120a 4 +32,256a 2 -663,552] 

(56) 


(55) 
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(a) Rectangular span loading. 

Figure 9. Power spectra of rolling moment due to horizontal gusts. 


A comparison of the values obtained for equations (55) 
and (56) and plots of the results obtained by the numerical- 
integration process indicated no difference, and none is shown 
in figure 9. 

It is significant to observe that very little variation exists 
in the power spectra of figure 9 for the four span loadings 
considered. However, as compared with the rolling moment 
due to vertical gusts (fig. 7), the rolling moment due to 
horizontal gusts is relatively small for small values of trim 
angle of attack. Although no exact expression for the ra t i o 
of the power spectra of the rolling moments due to u g and w g 
may be given without including /3' and y, it may be seen 
from figures 7 and 9 that, in general, 


0.2 


where a 0 is given in radians. 


(57) 


CALCULATIONS FOR SIDE GUSTS 


For the side gust considered, the correlation function of 
the rolling-moment coefficient as given by equation (32) 
becomes 


F, 


c,(Ut) 


U 2 


0 -^) 




(58) 



(b) Elliptic span loading. 
Figure 9. — Continued. 


equal to the variation of equation (41) with r, which is 
plotted in figure 4. 

The power spectrum of the rolling-moment coefficient as 
given by equation (33) with G(k’) given by equation (43) 
becomes 


$<?,(*') = 


Ci\v 2 L l +3(fc') 2 
tU 3 [1+(W 


(59) 


The variation of tire spectrum with frequency k'=~ is 
shown as the G(k') curve of figure 5. 

DISCUSSION 

The purpose of this section is to discuss the implications 
of the assumptions made in the analysis of this report, the 
reasons for making these assumptions, and the application 
of the results. 

ASSUMPTIONS CONCERNING THE NATURE OF TURBULENCE 

The turbulence was assumed to be homogeneous in order 
to make the problem stationary in the statistical sense and 
thus permit the use of the mathematical techniques developed 
for such problems. In a practical sense, turbulence can be 
homogeneous only in a limited body of air. The assumption 
thus implies that the dimension of this body of air along the 
flight path is large compared with the distance traversed 


and the variation of this function with t/r/L is, of course, 
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(c) Parabolic span loading. 
Figure 9. — Continued. 


in the reaction time of the airplane. In the case of loads 
studies this reaction time is of the order of the time to damp 
to one-half amplitude, but, in the case of motion studies, 
the reaction time may be much larger. Obviously, the 
greater the body of air, the greater the reliability with 
which the loads and motions can be predicted (in a statistical 
sense) for one run through it. In general, turbulence at 
very low altitudes, which may be influenced significantly by 
the configuration of the ground, and the turbulence in 
thunderstorms may not be sufficiently homogeneous for this 
type of analysis, but other types of turbulence are likely to 
be substantially homogeneous over sufficiently large distances. 

Isotropy was assumed in order to permit the required 
two-dimensional correlation functions to be expressed 
simply in terms of the one-dimensional correlation functions. 
For sufficiently short wavelengths all turbulence is isotropic 
(see ref. 7), but for long wavelengths it can be isotropic only 
if it is homogeneous (both in the plane of the flight path 
and perpendicular to it). (The condition of axisymmetry is 
less restrictive inasmuch as it does not specify the variation 
of the characteristics of the turbulence in the vertical 
direction.) In practical problems, if the turbulence may be 
assumed to be homogeneous, the conditions of isotropy are 
likely to be satisfied sufficiently to permit the use of the 
approach presented herein for all but very long wavelengths. 
The wavelength at which this approach ceases to be valid 



(d) Triangular span loading. 
Figure 9. — Concluded. 


depends on the size of the body of air under consideration, 
being longer for a large body. 

Taylor’s hypothesis implies that the variation in gust 
intensity that prevails along the flight path at any instant 
will remain substantially the same until the airplane has 
traversed the given body of air. The required correlation 
functions for atmospheric turbulence are thus in the nature of 
space correlation functions (rather than time correlation 
functions) and have been considered as such. The statistical 
characteristics of the turbulence are then independent of 
the speed at which it is traversed. Clearly, the validity of 
this hypothesis depends on the flying speed of the airplane 
and it would be expected that, at very low speeds, the 
hypothesis of Taylor becomes less valid and the results may 
be less accurate. On the basis of present knowledge, no 
definite lower limiting speed can be quoted. . The effect 
of finite flying speed on the gust correlation function can be 
expected to be most pronounced for large distances, where 
the correlation is weak. Thus, the effect on the various 
spectra is likely to be small and to occur at the longest 
wavelengths, where, as previously mentioned, the spectra 
are somewhat uncertain for other reasons as well. 

For practical purposes, the parameter L (the integral scale 
of turbulence) used herein is a largely fictitious quantity, 
inasmuch as it is, to a large extent, proportional to the values 
of the gust spectra for infinite wavelengths. In view of the 
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uncertainties in the values of the spectra at long wave- 
lengths and the fact that the spectra in this region pre- 
dominantly define the area under the integral, the parameter 
L has little physical significance. Therefore, at present, 
insufficient information is available to give an exact value 
for L to be used in connection with the numerical results 
calculated herein. However, on the basis of the measure- 
ments such as those of reference 13, a value of 1,000 to 2,000 
feet appears to be appropriate for the conditions of the 
referenced tests. It is desirable to obtain more information 
concerning the spectra of atmospheric turbulence under a 
wider range of conditions. More definite values could then 
be deduced by fitting measured results by means of an 
analytical expression of the type used here. This expression 
could be used as a means.of obtaining a value of L by extrap- 
olation of the measured results to infinite wavelengths (zero 
frequency) . 

ASSUMPTIONS CONCERNING THE AERODYNAMIC FORCES 

The fundamental assumption concerning the aerodynamic 
forces is that they vary linearly with gust intensity. This 
assumption implies that the ratio of the gust speed to the 
flying speed must always be fairly small; if the aerodynamic 
forces and moments tend to vary with gust intensity in a 
nonlinear manner, as the wing yawing moments do for all 
angles of attack and the other forces and moments do for 
high angles of attack, the ratio of gust intensity to flying 
speed must be very small — about 1/30 or less. However, 
as previously mentioned, the wing yawing moments due to 
gusts are likely to be quite small, so that some error in them 
due to slight deviations from linearity is not likely to affect 
appreciably the results of an analysis of the lateral motion. 
Hence, for an airplane flying at small angles of attack and 
at speeds of about 200 knots or more, in continuous turbu- 
lence, the assumption of linearity should lie valid; for flight 
in severe thunderstorms, it is not likely to be valid, and, for 
flight at high angles of attack, it is likely to be valid only for 
light turbulence. 

The rigidity of the wing, which was mentioned in the 
list of assumptions, enters only indirectly into the problem 
considered herein. The results obtained here are valid 
whether the wing is rigid or not. However, in the case of 
flexible wings (the term “flexible” being used to describe 
wings with deformations which give rise to appreciable 
aerodynamic forces), certain additional information is 
required. (See ref. 6 .) This information may take the form 
of span influence functions 7 (?/) modified by static aero- 
elastic effects, or may require certain cross-correlation 
functions or cross spectra between the gust forces and the 
dynamic forces, depending on the individual case. 

The assumption that the indicial-response influence func- 
tion h(t,y) can be written as a product of functions of time 
only and distance along the span only is based on the reason- 
ing of reference 6 . This reasoning, in turn, is based on the 
observation that, according to the available information for 
the lift distributions due to sinusoidal motions (and, hence, 
those due to indicial motion), the lift distribution tends to be 
substantially invariant with frequency (or time) except for an 
overall factor. Inasmuch as this information is confined to 


unswept wings, this assumption may not be valid for swept 
wings. 

APPLICATION OF THE RESULTS 

In this report the rolling moments and yawing moments 
have been calculated for a wing due to the u , v , and w com- 
ponents of turbulence. If the turbulence is isotropic, these 
components are statistically independent at a point. In any 
practical application, all three components are always present 
and the wing rolling and yawing moments due to the com- 
bined action of the three components must be known. In 
isotropic turbulence, the cross correlations between u and w 
and between v and w in the horizontal plane are zero, although 
u and v have a nonvanishing cross correlation. Thus, the 
moments due to v and w can be added directly, but, if hori- 
zontal-gust effects are to be taken into account, not only the 
moments due to u calculated herein but also the moments 
which arise from the cross correlation between u and v should 
be added to the others. However, there is reason to believe 
that the horizontal-gust effects on the lateral moments are 
generally very small, so that neglect of this cross-correlation 
effect is usually justified. 

The rolling and yawing moments due to u g and w g con- 
sidered herein are only those contributed by the wing but, 
inasmuch as the lateral moments contributed by the fuselage 
and tail as a consequence of these two components of gusts 
are generally very small, the results given here may, in 
general, be used to represent the lateral moments on a com- 
plete airplane due to these two gust components. 

Similarly, the rolling and yawing moments of a complete 
airplane due to the ^-component of gusts depend not only on 
the wing contribution considered here but also on the con- 
tribution of the vertical tail, which can be calculated in a 
straightforward manner. For instance, a method of calcu- 
lating the yawing moments and side force on a fuselage and 
vertical fin due to side gusts is found in reference 16. 

Although the contribution of the horizontal component of 
gusts to the lateral moments appears to be small compared 
with the other two components, it should be kept in mind 
that the effect of this component increases as the square of 
the trim angle of attack. (See eq. (57).) For conventional 
airplanes in the landing configuration and for vertically rising 
airplanes in the transitional stage, the effects of horizontal 
gusts may well be predominant in calculations of the forces, 
moments, and motions due to turbulence. 

CONCLUDING REMARKS 

The correlation functions and power spectra of the rolling 
and yawing moments 011 an airplane wing due to the three 
components of continuous random turbulence have been cal- 
culated. The rolling moments due to the longitudinal 
(horizontal) and normal (vertical) components depend on the 
spamvise distributions of instantaneous gust intensity, which 
were taken into account by using the inherent properties of 
symmetry of isotropic turbulence. The results consist of 
expressions for the correlation functions and spectra of the 
rolling moment in terms of the point correlation functions of 
the two components of turbulence. 
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Specific numerical calculations were made for a pair of 
correlation functions given by simple analytic expressions, 
which fit available experimental data very well. Calcula- 
tions were made for four lift distributions and the differences 
in the results calculated for these distributions were small. 
By comparison with the results calculated herein, the results 
of previous analyses for which it was assumed that random 
turbulence along the flight path and variations of turbulence 
across the span were linear have been shown to be valid only 
when the ratio of the span to the integral scale of turbulence 
(about 1,000 to 2,000 feet) is small. 

A comparison of the power spectra of the rolling moments 
due to horizontal gusts and those due to vertical gusts showed 
that the vertical gusts were predominant at small values of 


trim angle of attack (or trim lift coefficient) ; however, the 
relative effect due to horizontal gusts increased as a function 
of the square of the trim angle of attack. 

The rolling moment due to lateral (side) gusts, which is 
small, was expressed in terms of the instantaneous value of 
the gust at representative points on the wing, so that the 
effect of spanwise variation in gust intensity was ignored. 
The yawing moments were considered to be proportional to 
the rolling moments, the constants of proportionality being 
given by simple aerodynamic relations. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 6 , 1956. 


APPENDIX 


evaluation of the elliptic integral weighting function 

The evaluation of the integral weighting function T (77) involves the integral given by equation (21): 

r (rj) = J+ y(yf) y(yi*+ v) <ty* 

For the case of the elliptic distribution of the additional span loading factor, 

y(y*)='-^y\^-y * 2 


and the integral weighting function to be evaluated becomes 

r0?)=(^r) j i " ( y *) (y* -I- v) Vi — yi* 2 Vi — ( y*+v) 2 dyt * 


Under the substitution 


2y 1 »+iy __ 2 ? /,*+ij 
2-77 5 


the integral may be written as 

Inasmuch as 1 — z 2 = (1 — z)(l + 2), 

r(i7) 


(&r— 77) (5.2+77) -y/ 1 ^ (8x 77) 2 "y/ 1 ^ + 77) 2 dx 






dx 


8 2 — 77 

With the notation k=^-—=~ , > 


2+7? 2+7? 

r(7?)=+2 — 7?) 2 (2+t?) f [(2— 1 ?) V— 7? 2 ] V (1 —kx) (1 +x) (1 —a;) (1 +kx) 
TT J - 1 


dx 


=5 ( 2 “^) : 


,2++j; [(2-7?)V-7? 2 ] V(1 -k 2 x 2 ) (1-+) dx 
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where the integrand may be seen to be an even function of the variable x. Multiplying numerator and denominator by the 
radical and expanding yields 


T(v)=^(2- v y(2+v)S k\2- n y r -= f ' f ■ - - k>[(2+r,y+(2-r,y+r,*} P 

* k jo \(l—k 2 x 2 )(l—x 2 ) Jo 

[(2 - v y+v 2 +rk-] P . xHx =-y f 1 = = \ 

Jo V (1 —k' 2 x 2 ) (1 —x 2 ) Jo V(l-*V)(l-x 2 )/ 


x*dx 


V(l— k 2 x 2 )(l — x 2 ) 


The integrals may be recognized as elliptic integrals in powers of x 2n for which the closed-form solutions may be found in 
reference 1 7, for example. In terms of the standard elliptic integrals (in Jacobi's notation) of modulus k =-— —> 

2 + rj 


K(k)= f * , = 

Jo V(l— a:) (l—k 2 x 2 ) 


is defined as a complete elliptic integral of the first kind, and 



1 —k 2 x 2 

1— x 2 


dx 


is defined as a complete elliptic integral of the second kind. Tables of these integrals may be found in most mathematical 
handbooks as well as in reference 17. In terms of these integrals, the solution for the integral weighting function is found 
to be 


rfct)-“ (2- v y(2+r l ) [(8+3P+4 E)K(k) - (8+7F- +8^) E(k)]~ 

F[(2+^) 2 +(2-^) 2 +^] [{2+k * )K{k) - 2 (l+k>)Em + 

[K(k)-E(k)\-yK(Ic) j- 

512 

= 15^2 ( 2 +»)) [4i? (i? 3 — 3»7— 1 )K(k) + (4+9tj 2 — y)E(k)] 
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